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In dense neutron matter under the presence of a strong magnetic field, considered in the model 
with the Skyrme effective interaction, there are possible two types of spin ordered states. In one 
of them the majority of neutron spins are aligned opposite to magnetic field (thermodynamically 
preferable state), and in other one the majority of spins are aligned along the field (metastable 
state). The equation of state, incompressibility modulus and velocity of sound are determined in 
each case with the aim to find the peculiarities allowing to distinguish between two spin ordered 
phases. 
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I. INTRODUCTION 

Magnetars are strongly magnetized neutron stars [l| 
with emissions powered by the dissipation of magnetic en- 
ergy. Magnetars are thought to give the origin to the ex- 
tremely powerful short-duration 7-ray bursts 0, Q . The 
magnetic field strength at the surface of a magnetar is 
about of 10^^-10^^ gH,!!]- Such huge magnetic fields can 
be inferred from observations of magnetar periods and 
spin-down rates, or from hydrogen spectral lines. In the 
interior of a magnetar the magnetic field strength may be 
even larger, reaching values of about 10^* G [a0|- Under 
such circumstances, the issue of interest is the behavior 
of neutron star matter in a strong magnetic field ■ 

In the recent study 0, neutron star matter was ap- 
proximated by pure neutron matter in a model with the 
effective nuclear forces. It was shown that the behavior 
of spin polarization of neutron matter in the high den- 
sity region in a strong magnetic field crucially depends on 
whether neutron matter develops a spontaneous spin po- 
larization (in the absence of a magnetic field) at several 
times nuclear matter saturation density, or the appear- 
ance of a spontaneous polarization is not allowed at the 
relevant densities (or delayed to much higher densities). 
The first case is usual for the Skyrme forces [lol - [2l| , while 
the second one is characteristic for the realistic nucleon- 
nucleon (NN) interaction [2^ - [28j . In the former case, a 
ferromagnetic transition to a totally spin polarized state 
occurs while in the latter case a ferromagnetic transition 
is excluded at all relevant densities and the spin polariza- 
tion remains quite low even in the high density region. 
If a spontaneous ferromagnetic transition is allowed, it 
was shown in the subsequent model consideration with 
the Skyrme effective forces [l^ that the self-consistent 
equations for the spin polarization parameter at nonzero 
magnetic field have not only solutions corresponding to 
negative spin polarization (with the majority of neutron 
spins oriented opposite to the direction of the magnetic 



field) but, because of the strong spin-dependent medium 
correlations in the high-density region, also the solutions 
with positive spin polarization. In the last case, the for- 
mation of a metastable state with the majority of neutron 
spins oriented along the magnetic field is possible in the 
high-density interior of a neutron star. 

In the present study, we provide the zero-temperature 
calculations of the equation of state (EoS), incompress- 
ibility modulus and sound velocity for neutron matter 
in a strong magnetic field with the aim to find the pecu- 
liarities allowing to distinguish between two possible spin 
ordered states - the stable one with negative spin polar- 
ization and the metastable one with positive spin polar- 
ization. It will be shown that in the thermodynamically 
stable state the incompressibility modulus and the speed 
of sound arc characterized by the appearance of the well- 
defined maximum just around the density at which the 
ferromagnetic (FM) phase transition sets in. Contrarily 
to that, such features are missing in the metastable state. 
Besides, all calculated quantities behave differently under 
changing magnetic field in stable and metastable states. 

At this point, it is worthy to note that we consider ther- 
modynamic properties of spin polarized states in neutron 
matter in a strong magnetic field up to the high den- 
sity region relevant for astrophysics. Nevertheless, we 
take into account the nucleon degrees of freedom only, 
although other degrees of freedom, such as pions, hyper- 
ons, kaons, or quarks could be important at such high 
densities. 



II. BASIC EQUATIONS 

Here we only outline the basic equations necessary for 
further calculations, and a more detailed description of 
a Fermi-liquid approach to neutron matter in a strong 
magnetic field can be found in our earlier work [2^ . The 
normal (nonsuperfluid) states of neutron matter are de- 
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scribed by the normal distribution function of neutrons 
/kik2 — Trgia^^aKi, where k, = (p, cr), p is momentum, 
a is the projection of spin on the third axis, and g is 
the density matrix of the system [l^, [l^ [2l| . Further 
it will be assumed that the third axis is directed along 
the external magnetic field H. The self-consistent matrix 
equation for determining the distribution function / fol- 
lows from the minimum condition of the thermodynamic 
potential [s^l and is 



/ = {exp(ro£ + + 1} ' = {exp(yoO + 1}" 



(1) 



Here the single particle energy e and the quantity Y/^ 
arc matrices in the space of k variables, with = 
5^4^KiK2' ^0 = ^/T, and I4 = — Ato/T being the Lagrange 
multipliers, /ip being the chemical potential of neutrons, 
and T the temperature. Given the possibility for align- 
ment of neutron spins along or oppositely to the magnetic 
field H, the normal distribution function of neutrons and 
single particle energy can be expanded in the Pauli ma- 
trices fji in spin space 



/(P) = /o(p)o-o + /3(P)<73, 

e(p) = £o(p)o-o +e3(p)cr3. 



(2) 



Using Eqs. ([T} and one can express evidently the 
distribution functions /o, /a in terms of the quantities e: 



/o = \{n{i^+)+n{LU-)}. 
h = \{n{i^+) - n{Lu^)}. 



Here n{uj) = {exp(yow) + 1} ^ and 

^± = Co ±6, 
Co = £0 - A^o, = £3- 



(3) 



(4) 



As follows from the structure of the distribution func- 
tions /, the quantities io± play the role of the quasipar- 
ticle spectrum and correspond to neutrons with spin up 
and spin down. The distribution functions / should sat- 
isfy the normalization conditions 
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y'^hiv) = Q\- Qi = A(?. 



(5) 
(6) 



Here g = + Qi is the total density of neutron matter, 
g^ and g\^ are the neutron number densities with spin up 
and spin down, respectively. The quantity IS.Q may be 
regarded as the neutron spin order parameter. The spin 
ordering in neutron matter can also be characterized by 
the neutron spin polarization parameter 



H = 



The spin order parameter determines the magnetization 
of the system M = n^Ag, fi^ being the neutron mag- 
netic moment. The magnetization may contribute to 
the internal magnetic field B = H + AttM. However, 
we will assume, analogously to Refs. 0, 0; that the 
contribution of the magnetization to the magnetic field 
B remains small for all relevant densities and magnetic 
field strengths, and, hence, B H . This assumption 
holds true due to the tiny value of the neutron mag- 
netic moment /i„ = -1.9130427(5)Aijv ~ -6.031 • IQ-^^ 
MeV/G [m (/iAT being the nuclear magneton) and is con- 
firmed numerically in a subsequent integration of the self- 
consistent equations. 

In order to get the self-consistent equations for the 
components of the single particle energy, one has to set 
the energy functional of the system. In view of the above 
approximation, it reads [l9| 

E{f) = Eoif, H) + E.mU) + Ef,M, (7) 

i?o(/,i?) = 2^£o(p)/o(p)-2/i„i7^/3(p), 
p p 

E^ntif) = 5^{eo(p)/o(p) + £3(P)/3(P)}, 

p 

i/2 



E field = "^V, 
OTT 



where 



£o(p) = ^Ef^o(k)/o(q), k = 



p q 

2 ' 



£3(p)-^Et^r(k)/3(q)- 



(8) 



(9) 



Here e;o(p) = is the free single particle spectrum, 
mo is the bare mass of a neutron, [/^(k), t/f (k) are the 
normal Fermi liquid (FL) amplitudes, and eo,e3 are the 
FL corrections to the free single particle spectrum. Note 
that the field contribution Efi^u, being the energy of the 
magnetic field in the absence of matter, leads only to the 
constant shift of the total energy and, by this reason, can 
be omitted. Using Eq. 0, one can get the self-consistent 
equations in the form [l9| 



Co(p) = £o(p) +£o(p) - Mo, 
C3(P) = -/i„-ff + e3(p)- 



(10) 
(11) 



To obtain numerical results, we utilize the effective 
Skyrme interaction [32[. Expressions for the normal FL 
amplitudes in terms of the Skyrme force parameters were 
written in Refs. [13, [l!]. Thus, using expressions ([3]) for 
the distribution functions /, we obtain the self-consistent 
equations ([TUl) . (|lip for the components of the single- 
particle energy Co(p) and CsIp), which should be solved 
jointly with the normalization conditions ([5|), (jB]). Fur- 
ther we do not take into account the effective tensor 
forces, which lead to coupling of the momentum and spin 
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degrees of freedom, and, correspondingly, to anisotropy 
in the momentum dependence of FL amplitudes with re- 
spect to the spin quantization axis. 

If the self-consistent equations have a few branches of 
the solutions, it is necessary to compare the correspond- 
ing energies (at zero temperature) in order to decide 
which solution is thermodynamically preferable. The en- 
ergy per neutron, E/A, can be directly calculated from 
Eq. The equation of state (EoS) of neutron matter 
in a strong magnetic field then can be obtained from the 
equation 



P 



dg 



(12) 



where e = g{mc^ + E/A) is the energy density, which 
includes also the rest energy term. The incompressibility 
modulus, K = 9^, according to Eq. (|T2|) . reads 



,d^E/A) P 
K = , ' + 18-. 



Q 



(13) 
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FIG. 1: Pressure vs. density for the branches Hi (stable) and 
Ilg (metastable) of spin polarization in neutron matter with 
the Skyrme SLy7 interaction aX U = 10^* G. 



The speed of sound, = cy can be related to the 
incompressibility modulus by the equation 



K 



+ 7) 



(14) 



III. 



EOS OF DENSE NEUTRON MATTER IN A 
STRONG MAGNETIC FIELD 



The self-consistent equations were analyzed at zero 
temperature in Ref. [2^ for the magnetic field strengths 
up to Hmax ~ 10^^ G, allowed by a scalar virial the- 
orem |34| . in the model consideration with SLy4 and 
SLy7 Skyrme effective forces [s^. It was shown that 
a thermodynamically stable branch of solutions for the 
spin polarization parameter as a function of density cor- 
responds to negative spin polarization when the majority 
of neutron spins are oriented opposite to the direction of 
the magnetic field. Besides, beginning from some thresh- 
old density gth ~ 4giQ, being slightly dependent on the 
magnetic field strength, the state with positive spin po- 
larization can also be realized as a metastable state in 
neutron matter (cf. branches Hi and li^ in Fig. 2 of 
Ref. [2^ ) . This conclusion was based on the comparison 
of the free energies of two states which turn out to be 
very close to each other [2^, [s^. However, as it will be 
shown later, additional constraints, such as, e.g., stabil- 
ity of the system with respect to the density fluctuations, 
will define more accurately the density range admissible 
for the state with positive spin polarization. 

In this work, the previous study [29j will be extended 
by calculating the EoS of dense neutron matter in a 
strong magnetic field for various branches of solutions of 
the self-consistent equations. Each possible state should 



match the constraint > for allowable densities and 
magnetic field strengths being the condition of the me- 
chanical stability of the system. Besides, in the high- 
density region, the velocity of sound should not exceed 
the speed of light in the vacuum, Vs < c. Note that fur- 
ther the contribution of the magnetic field pressure to 
the total pressure will be omitted because in the mag- 
netic fields up to 10^® G the magnetic field pressure is 
still small compared to the matter pressure in the high- 
density region of interest. 

First, we present the results of determining the zero- 
temperature EoS of neutron matter in a strong mag- 
netic field at the density region where both stable and 
metastable spin ordered states can be realized. Because 
the results of calculations with SLy4 and SLyT Skyrme 
forces are very close, here we present the obtained de- 
pendences only for the SLy7 Skyrme interaction. Fig. 1 
shows the pressure of neutron matter as a function of den- 
sity for two branches of spin polarization, stable Hi and 
metastable Ha, corresponding to negative and positive 
polarizations, respectively (the branch 112 with positive 
spin polarization considered in Ref. p9j has the consider- 
ably larger energy per neutron as compared to the previ- 
ous ones). For the branch Hi, the pressure is the increas- 
ing function of the density for all relevant densities, and, 
hence, the incompressibility coefficient is always positive. 
However, for the branch H3, beginning from the thresh- 
old density gth ~ 3.92£io up to the density gc ~ 4.85po 
(at H = 10^^ G), the pressure decreases with the density. 
Hence, in this density range the incompressibility coeffi- 
cient is negative and the metastable state characterized 
by the branch H3 of positive spin polarization cannot 
appear at these densities. However, beyond the critical 
density g^ the metastable state with positive spin po- 
larization is allowed by the criterion K > 0. Note that 
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FIG. 2: Pressure vs. magnetic field strengtli for tfie brandies 
III and Us of spin polarization in neutron matter with the 
Skyrme SLy7 interaction at g — 5go. 



FIG. 3: Same as in Fig. 1, but for the incompressibility mod- 
ulus of neutron matter in a strong magnetic field. 



the EoS for the mctastablc state of neutron matter in a 
strong magnetic field is stiffer than that for the thermo- 
dynamically equilibrium state. 

Fig. 2 shows the pressure of neutron matter as a func- 
tion of the magnetic field strength for the branches Hi 
and Ha of spin polarization at g = 5go- It is seen that 
the dependence of the EoS for stable and metastable 
branches is different: for the branch Hi the EoS becomes 
softer with the magnetic field while for the branch Ha 
stiffer. These calculations show that the impact of the 
magnetic field on the EoS remains small up to the field 
strengths of about 10^^ G. 

Fig. 3 shows the zero-temperature incompressibility 
modulus of neutron matter in a strong magnetic field 
as a function of density for the branches Hi and lis of 
spin polarization. For the branch 1X3, the incompress- 
ibility modulus monotonously increases with the density 
and changes sign from negative to positive at the critical 
density Qc, marking the stability range with respect to 
density fluctuations at densities beyond Qc- As a con- 
sequence, if the metastable state with positive spin po- 
larization can be realized in the high-density region of 
neutron matter in a strong magnetic field, under decreas- 
ing density (going from the interior to the outer regions 
of a magnetar) it changes at the critical density Qc to a 
thermodynamically stable state with negative spin polar- 
ization. 

For the branch Hi, the behavior of incompressibility 
modulus is nonmonotone. The most important pecu- 
liarity is that just around the density {gtr ~ 3.16^0 at 
H = 10^^ G) at which the magnitude of the spin po- 
larization parameter for the branch IIi begins rapidly 
to increase (cf. Fig. 2 of Ref. H^), the increasing be- 
havior of the incompressibility modulus with the density 
changes on the decreasing one. Because the density gtr 
can be regarded as the density at which a ferromagnetic 



state sets in, this qualitative feature in the behavior of 
the incompressibility modulus can be used as the char- 
acteristic of the density-driven FM phase transition in 
neutron matter possessing equilibrium spin polarization. 
The incompressibility modulus decreases till the density 
about 4po at which the spin polarization parameter is 
well developed and gets about two third of its strength. 
Then there is the plateau in the density dependence of 
incompressibility modulus till the density about 5ga be- 
yond which the incompressibility modulus begins grad- 
ually to increase. Note that the noticeable decrease of 
the incompressibility modulus around the density of the 
transition to the ferromagnetic state was mentioned also 
in Ref. [13] , although the total incompressibility modulus 
was not explicitly shown there, but only that for the spin- 
up and spin-down neutron components in the state with 
equilibrium spin polarization. Besides, in Refs. jol. [37j. 
there were no any calculations related to the metastable 
branch of spin polarization because this branch itself was 
missed in these studies. 

Fig. 4 shows the incompressibility modulus of neutron 
matter as a function of the magnetic field strength for 
the branches Hi and Ha of spin polarization at the den- 
sity g = 5go. For the branch Hi, the incompressibility 
modulus increases with the magnetic field strength while 
for the branch 1X3 it decreases. It follows from these cal- 
culations that the impact of the magnetic field on the 
incompressibility modulus remains mild up to the field 
strengths of about 10^'' G. 

Fig. 5 shows the sound velocity in neutron matter un- 
der the presence of a strong magnetic field as a func- 
tion of density for the branches Hi and II3 of spin po- 
larization. For the branch 1X3, Eq. ([T?|) automatically 
guarantees the fulfillment of the condition K > (for 
all relevant densities E/A > 0). It is seen that for 
both branches at the relevant densities the superluminous 
regime doesn't occur. While for the branch Ha the sound 
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FIG. 4: Incompressibility modulus vs. magnetic field strength 

for the branches Hi and of spin polarization in neutron pj^ 5. ^^^^ ^ p^g ^ut for the sound velocity in neu- 
matter with the Skyrme SLy7 mteraction at g = 5^0. tron matter under the presence of a strong magnetic field. 



velocity monotonously increases with the density, for the 
branch Hi it has non-monotone behavior. In fact, near 
the transition density gtr the sound velocity in a ther- 
modynamically stable state has a clear peak structure 
considerably decreasing at the densities where the ferro- 
magnetic phase sets in. This feature, together with the 
presence of the maximum in the density dependence of 
the incompressibility modulus, can be used for the iden- 
tification of the density-driven FM phase transition in 
neutron matter possessing equilibrium spin polarization. 
On the other hand, the incompressibility modulus and 
the speed of sound monotonously increase with density 
in the metastable state with positive spin polarization, 
and hence, these features can be used for distinguish- 
ing between the thermodynamically stable (negative spin 
polarization) and metastable (positive spin polarization) 
states in neutron matter under the presence of a strong 
magnetic field. 

Fig. 6 shows the sound velocity in neutron matter as a 
function of the magnetic field strength for the branches 
Hi and 113 of spin polarization at the density g = 5go- 
For the branch Hi, the sound velocity increases with the 
magnetic field strength while for the branch II3 it de- 
creases. These trends are quite similar to those in the 
behavior of the incompressibility modulus K{H) for the 
branches Hi and 



IV. CONCLUSIONS 

In dense neutron matter under the presence of a strong 
magnetic field, considered in the model with the Skyrme 
effective interaction, there are possible two types of spin 
ordered states: the one with the majority of neutron spins 
aligned opposite to magnetic field (thermodynamically 
preferable state), and the other one with the majority 
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FIG. 6: Same as in Fig. 4 but for the sound velocity in neutron 
matter under the presence of a strong magnetic field. 

of spins aligned along the field (metastable state). The 
equation of state, incompressibility modulus and veloc- 
ity of sound have been determined in each case for SLy7 
Skyrme force with the aim to find the peculiarities allow- 
ing to distinguish between two spin ordered phases. 

For the stable state with the branch Hi of negative spin 
polarization, the EoS is softer than that for metastable 
state with the branch II3 of positive spin polarization. 
The condition of the positiveness of the incompressibility 
modulus, K > 0, is satisfied for all relevant densities and 
magnetic field strengths for the stable branch Hi. How- 
ever, for the branch although formally the solutions 
of the self-consistent equations exist at densities larger 
than some threshold one, gth^ the condition K > is 
satisfied only at the densities larger than the critical one, 
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Qc (e.g., for H = G, Qth « 3.92^0 and Qc ~ A.Sbgn). 
As a consequence, if the metastablc state with positive 
spin polarization can be reahzed in the high-density re- 
gion of neutron matter in a strong magnetic field, under 
decreasing density (going from the interior to the outer 
regions of a magnctar) it changes at the critical density Qc 
to a thermodynamically stable state with negative spin 
polarization. 

For the thermodynamically stable branch Hi, the in- 
compressibility modulus and the speed of sound are char- 
acterized by the appearance of the well-defined maximum 
just around the density at which the ferromagnetic phase 
sets in. The last qualitative features can be used for the 
identification of the density-driven FM phase transition 
in neutron matter, possessing equilibrium spin polariza- 
tion, under the presence of a strong magnetic field. Con- 
trarily to the previous case, for the branch II3 of positive 
spin polarization, the incompressibility modulus and the 
speed of sound monotonously increase with density that 
can be used to distinguish between two different spin or- 
dered phases. 



The dependence of all calculated quantities on the 
magnetic field strength H turns out to be different for two 
spin ordered phases. For the thermodynamically stable 
branch Hi , the incompressibility modulus and sound ve- 
locity increase with H while the pressure decreases. The 
exactly opposite tendency has been found for the branch 
II3 of positive spin polarization that also allows one to 
differentiate between spin polarized states with opposite 
polarizations in neutron matter under the presence of a 
strong magnetic field. 

As yet one problem for consideration, it would be in- 
teresting to study the role of finite temperature effects on 
the EoS, incompressibility modulus and speed of sound in 
dense neutron matter in a strong magnetic field. As has 
been shown already, these effects can lead to a number 
of nontrivial features such as, e.g., unusual behavior of 
the entropy in various spin ordered systems [s^ Issl - lloj . 
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MEST and NRF. 
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